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1. Introduction 

Ionization is one of the most interesting physical 
processes arising from the interaction of atoms 
or molecules with electromagnetic fields. Ionized 
electrons are characterized by their ionization spectra 
that contain information about the electronic structure 
of the ionized atoms or molecules. For this reason, 
ionization processes have been used as a strong tool 
for the investigation of electronic structures (TJ [5] . The 
ionization spectra are strongly sensitive to individual 
bound excited electronic states through which the 
ionization can efficiently proceed (for an extended list 
of references, see e.g. [a SI 016]). 

The mechanism of auto-ionization occurring with 
the contribution of a bound excited state has been 
extensively studied for He atoms containing only two 
electrons. Quantum mechanical calculations have 
shown that there is a large probability of exciting 
both electrons together [7] due to the strong electron- 
electron correlations present in the electronic ground 
state [H|. Energy of one of the excited electrons given 
to the other excited electron then allows ionization of 
this electron. The first electron returns to its ground 
state. This process modelled in [a [10 [m with the 
inclusion of an additional probe field has recently been 
observed using attosecond pump-probe spectroscopic 
techniques [laiia- 

The presence of an auto-ionizing state in atoms 
manifests itself in the ionization spectra by the so- 
called Fano zero [14]. It denotes an electron energy 
missing in the ionization spectrum. This occurs due 
to the completely destructive interference between 
different paths of ionization, including the interference 
between direct and indirect (through an auto-ionizing 
state) ionization processes [15]. This interference also 
results in sharp peaks in the ionization spectra caused 
by intense laser pulses [Min]- As predicted by Fano 
|14j . n zeroes will be present for n auto-ionizing levels 
for an atom as long as only a single ionization channel 
is available. They are replaced by local minima when 
more ionization channels are present, and they may 
even disappear in molecules because of the presence of 
nuclear motion |18j . 

The presence of auto-ionizing states influences 
considerably many physical effects. For example, auto¬ 
ionization systems have been investigated as media 
exhibiting electromagnetically-induced transparency 
[laiioiiiT]. They have also been applied for slowing 


down the light propagating through the medium with 
auto-ionizing states [20]. Moreover, special attention 
has been devoted to the problem of near-threshold 
ionization [22] including quantum anti-Zeno effect |23j . 
The role of elastic collisions and finite laser spectral 
width has been analyzed for such systems in |24j . 
Auto-ionization has even been considered for weak 
quantum fields [25], in particular, for squeezed light 
pSj . Interestingly, Fano profiles and Fano resonances 
can appear not only in typical atomic/optical models. 
They have been found, for instance, in plasmonic 
systems m [28] . various nanoscale structures [29] 
including quantum dots [30] or in a broad range of 
superconducting systems [SI1I31I33]- We note that 
states of the continuum and bound states participate 
together in the effect called the Feshbach resonance 
[34] . In this effect, a quasi-bound molecular state is 
formed from the state of two free atoms. This results 
in considerable enhancement of inelastic collision 
processes [55] . 

Ionization can also be observed due to the 
electronic Coulomb interaction between two excited 
electrons being at two neighbor atoms. In this case, 
energy of one excited electron is given to the excited 
electron on the neighbor atom at which ionization 
occurs |36j . If two excited electrons emerge from 
the interaction with the strong electromagnetic field, 
we have two-center resonant photo-ionization [371138] . 
Also two-center dielectronic recombination [39l [38] 
which is just the inverse process in which a free electron 
is captured by a positively charged system of two 
neighbor atoms has been suggested. These effects are 
important in atomic and molecular clusters in which 
the exponential character of these processes can be 
modified by additional interactions, e.g. by capturing 
the ionized electron [40] . 

In an intense external field, double ionization may 
occur leaving two free ionized electrons [7] . It has been 
extensively studied in isolated atoms or molecules [41] . 
The dynamics of double ionization is predominantly 
governed by the strong Coulomb interaction among the 
electrons and the core of an atom. As a consequence, 
strong correlations in momenta of two ionized electrons 
occur [42l [43]. Considering He atoms, rotational 
and bending modes of the two-excited-electron motion 
preceded the ionization have been identified [23] . 
Double ionization of Ne atoms has been experimentally 
studied in [45] . 

Here, we consider another process in which 
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two electrons are ionized in such a way that their 
energies are highly correlated. Contrary to the double 
ionization model of a single atom, we consider a 
system of two neighbor atoms or molecules under the 
interaction with a strong cw external field. Each 
atom or molecule provides one ionized electron. We 
assume here that each atom or molecule can be ionized 
through an auto-ionizing level or by the process of 
direct ionization. We show that, as a result of 
the dipole-dipole interaction between two atoms or 
molecules, the state of two ionized electrons can be 
strongly entangled and both electrons can exhibit 
strong spectral correlations. 

The model considered here is analyzed using the 
algebraic method applied already by Fano [Tl]. Our 
model represents a generalization of the previously 
developed models in which ionization spectra and 
correlations between an ionized electron and a bound 
electron found at a neighbor atom have been studied. 
It has been shown that the so-called dynamical 
Fano zeroes occurring once per the Rabi period 
can be observed in conditional ionization spectra, 
independently on the presence [46] or absence [47] of 
an auto-ionizing level. Moreover, for the previously 
discussed systems, the dipole-dipole interaction leads 
to the generation of entangled states encompassing 
both the ionized and the bound electrons |48l [49] . The 
observed spectral correlations have been conveniently 
quantified by the density of quadratic negativity [bOj 
defined in terms of partially transposed statistical 
operators [511 [52]. Moreover, it has been shown that 
the dipole-dipole interaction does not suppress the 
occurrence of the Fano zero in general, but its existence 
is restricted to only rather special conditions [53l [54] . 
Also a quantized field has been considered for the two- 
atom auto-ionizing model showing that many sharp 
peaks in the ionization spectra originating in discrete 
field energies can occur [49]. 

It will be shown that the exact Fano zeroes typical 
for ionization spectra of individual atoms cannot occur 
in the analyzed model of two auto-ionization systems. 
On the other hand, many spectral features found 
in the Fano model of an auto-ionization system |14| 
are present in the analyzed system including sharp 
peaks and dips. The applied algebraic approach 
provides analytical formulas that allow to get a deep 
insight into the behavior of two interacting auto¬ 
ionization subsystems analyzing spectral correlations 
and entanglement [55] of the state of two ionized 
electrons. Similarly as in [48] the density of quadratic 
negativity is applied to quantify quantum spectral 
correlations. Although the model is discussed in 
general using the values of parameters emphasizing its 
main features, the values of parameters suitable for 
molecular condensates (including molecular crystals) 



Figure 1. Scheme of the mutually interacting auto-ionization 
systems a and b. State |l)j denotes an excited bound state of 
atom j with energy Ej whereas state \Ej) with energy Ej lies 
inside the continuum of atom j, j = a, 6. Symbols fij and fij 
stand for the dipole moments between the ground state |0)j and 
the corresponding excited states, ctL stands for the pumping 
amplitude, Vj describes the Coulomb configurational coupling 
between states |l)j and \Ej). Symbols Ja ^-nd Jf, refer to 
the dipole-dipole interactions between the excited/ionized states 
at atoms a and b indicated in the scheme. Double arrows mean 
that two electrons at atoms a and b participate in the interaction 
leading to energy transfer. 

[55] are also considered. 

The paper is organized as follows. Hamiltonian 
of two coupled auto-ionization systems in a stationary 
optical field is given in Sec. 2, together with the 
corresponding dynamical equations and their solution. 
Quadratic negativity and its density as quantifiers 
of entanglement between two ionized electrons is 
introduced in Sec. 3. Long-time photoelectron 
ionization spectra for the cases of comparable direct 
and indirect ionization paths (Subsec. 4.1) and 
dominating indirect ionization path (Subsec. 4.2) are 
analyzed in Sec. 4. Entanglement in the long-time 
photoelectron spectra is discussed in Sec. 5. Sec. 6 
is devoted to the process of ionization in molecular 
condensates. Sec. 7 brings conclusions. Analytical 
solution of the model without the mutual interaction is 
provided in Appendix A. Appendix B is devoted to the 
competition of the dipole-dipole interactions between 
discrete auto-ionizing levels and with the continua. 

2. Hamiltonian, dynamical equations and their 
solution 

We consider two auto-ionization systems describing, 
e.g., two atoms or molecules that mutually interact 
by the dipole-dipole interaction (for the scheme, see 
Fig.dj. Both the ionization systems denoted as a and 
b are assumed to have one auto-ionizing level. They 
interact with a stationary optical field with amplitude 
ul and frequency through the corresponding dipole 
moments. Hamiltonian iJj describing auto-ionization 
system j, j = a, 6, can be written as (?i = 1 is assumed, 

my 

H; =H° + H^, j = a,h, (1) 
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=E°\1),,{1\ + J dE,E,\E,){E,\ 

+ J dE,[V,\E,),{l\+R.c.], 

^^ljaL exp{-iELt)\l)jj{0\ 

+ / dEj fljaL exp{—iELt)\Ej) j{0\+li.c. 


H^ = 


( 2 ) 

(3) 


In Eqs. (P-El), Hamiltonian ijP describes the 
inner structure of auto-ionization system j whereas 
Hamiltonian arises from the interaction of the 
system with the optical field. Symbol Ej denotes 
the excitation energy from the ground state |0)j 
into the excited bound state |l)j of atom j. The 
continuum of atom j is formed by states \Ej) having 
energies Ej. The Coulomb configuration interaction 
between states |l)j and \Ej) inside the continuum is 
described by coupling constant Vj. Dipole moments 
^j and flj characterize an optical excitation of the 
corresponding states. Symbol H.c. stands for the 
Hermitian conjugated term. 

The dipole-dipole interaction between atoms a and 
b leads to energy transfer [56]. It can be described by 
the following Hamiltonian i/trans^ 


Htrpiris — 


Ja&|0)a|l)&h(0|a{I| 

+ J dEaJa\Ea) |0)6h(l|a{0| 

+ [ dEbJb\0}a\Eb)b{0\a{l\+R.C. 


(4) 


In this interaction, one electron loses its energy and 
returns into its ground state, whereas the other 
electron takes this energy and moves into its own 
excited/ionized state. Constant Jab describes the 
interaction between excited states |l)a and |1){, whereas 
constant Ja [Jb] characterizes the interaction between 
excited state |l)b [|l)a] and ionized states \Ea) [|Eb)]. 
We assume that the electrons in ionized states do not 
mutually interact. 

A general quantum state of two electrons at atoms 
a and b written in the frame rotating with frequency 
El can be expressed as: 


\tp){t) = Coo(t)|0)a|0)f, -f Cio(f)|l)a|0)b 

+ Coi(t)|0)a|I)b + Cll(t)|I)a|I)& 

+ J dEa [da,o{Ea, t)|£'a)|0)6 + da.l{Ea, t)|£'a)|l)b] 
+ J dEb[dbfl{Eb,t)\0)a\Eb) + db^i{Eb,t)\l)a\Eb)] 

+ J dEa J dEbd{Ea,Eb,t)\Ea)\Eb). (5) 


Time-dependent coefficients cqo, ciq, cqi, ch, 
^a,o(A'a), da.i(^Ea\ db.o{^Eb)j and di^Eaj Eb^ 

characterize the state \tp) at an arbitrary time 


t. They can be conveniently grouped into the 
vectors c^{t) = [coo(t), cio(t), coi(t), cii(t)] and 

dj{Ej,t) = [dj^o{Ej,t),dj^i{Ej,t)] for j = a,b. 
Symbol T denotes transposition. Normalization 
of the state |V^)(t) in Eq. (O provides the re¬ 
lation ct(t)c(t) -h 'Ej=a,bl dEjd]{Ej,t)dj{Ej,t) + 
f dEaf dEb\d(Ea,Eb,t)\^ = l. 

The Schrodinger equation with the overall Hamil¬ 
tonian Ha -f Hb -f I?trans provides the following dynam¬ 
ical equations for the coefficients characterizing state 
written in Eq. ([5)): 


c{t) 


A 

/ dEoBa 

da{Ea,t) 


Bi 

{Ea — El)1 -f Kb 

db{Eb,t) 


Bl 

0 

_ d{Ea,Eb,t) 


0 



/ dEbBb 

0 


c{t) 

0 

/ dEblb 


da{Ea,t) 

{Eb — El)1 + Kq 

f dEala 


db{Eb, t) 

li 

Ea -f Eb — 2 El 


_ d{Ea,Eb,t) 


( 6 ) 


In Eq. (|6|), symbol 1 (0) stands for the unity (zero) 
matrix of appropriate dimension(s). 

Matrices A, B^, Bf,, K^, Kf,, and lb introduced 
in Eq. (0) are defined as follows: 


A = 


0 

f^>*L 

* ^ 
d-b<^L 

0 



Jib 

* sk 

f^bO^L 

d'bOLL 

Jab 

AEO 


0 

^bOCL 

f^a^L 

AEl + AEl 


K, = 



' ~ 5k 5k 

t^a^L 

1 

O 


“ ~ 5k 5k 

f^bO^L 

1- 

o 

Ba = 

V* 

* a 

T* 

0 

~ 5k 5k 

, Bb = 

j: 

v: 

~ 5k 5k 

0 


0 

V* 

O. 


0 

-1 

* 

A" 


0 

^ijaL 


AE^ 


T — 

■ ■ 

1 

H 

A- 


j=a,b-, (7) 

AE° = E^-EL,j = a,b. 

The set of equations ([HI) can be conveniently solved 
by the Laplace transform defined for an arbitrary 
function f{t) as /(e) = dtf{t) exp{iet). If we 
assume that only the discrete levels of atoms a and 
b are initially populated, we arrive at the following set 
of equations: 


[el-A]c(e)- j dEj'Rjdj{Ej,e)=ic{Q), (8) 

i=a,bd 

[(e — Ea + El ) 1 — K{,] (Eq , e) 

- B))c(e) - J dEblbd{Ea, Eb,e) = 0, (9) 

[(£ — Eb + El)1 — Kq] db{Eb, e) 

-Blc{e) - J dEalad{Ea,Eb,e) = 0, (10) 
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(e — Ea — Eh + 2E]^)d{Ea, Eh, s) 

-lldaiE,,e)-lidh{Eh,e) = 0. ( 11 ) 

Vector c(0) occurring in Eq. dH) describes the initial 
condition. If the electrons at both atoms are initially 
in their ground states, c^(0) = [1, 0,0, 0]. 

Equations (EHni) can be solved as follows. First, 
using Eq. (fTTl) we express function d{Ea, Eh,e) in the 
following form: 

Ibda(£'ai e) + s) 


d{Ea, Eh,e) — 


( 12 ) 


£ — Ea — Eh + 2E]^ 

Snbstitntion of the expression (fT^ into Eqs. and 
m and subseqnent integration in these equations 
results in the following equations: 

[(& — Ea + El )1 + Lb] da{Ea, s) 
lhllMEb,e) 


= Bt 


C £ 


(13) 


= b1c(£). (14) 


- / dEh- 

J e-Ea-Eh + 2EL 

[(e — Eh + El) 1 + Jja] dh{Eh, e) 

- [dE^ I^lId.(K,£) 

J ^e-Ea-Eh + 2EL 

The newly defined 2x2 matrices and Lf, are given 
as 

hj = —Kj + iTrljlj, j = a, b. (15) 

In the next step, we neglect the terms in Eqs. m 
and da containing the integrations over Ea and Eh- 
As shown in Appendix A these terms equal zero for 
independent atoms a and b. Weakness of the dipole- 
dipole interaction compared to the Coulomb one then 
justifies this approximation that leaves us the following 
decoupled equations: 

[{e-Ea+EL)l + Lh]daiEa,e)=Blc{e), (16) 

[{£-Eh + EL)l + ha]dhiEh,e) =bJc(£). (17) 

These equations can be solved using the decomposition 
of matrices and u m, 


Ej = j = 


(18) 


k=l,2 


in which A'(. are eigenvalues of the matrix hj and 
matrices Lj and L 2 fulfil the relation Lj -f L 2 = 1. 
Formulas for eigenvalues Xj. are derived as follows: 


•^ 1,2 = «i ± V af - j = a,b; (19) 

ai = [z7r{\f,,aLf + \V,f)-AE^] /2, 
a^2 = -iTrlfijaLl^AEj - 

-I- ITT {fl*fj.j\aL\^Vj -I- c.c.) , 

where c.c. denotes the complex conjugated term. Using 
eigenvalues X^,, the matrices L-^, are expressed as: 


Substitution of the decomposition (TTSl) into Eq. dT51) 
allows to obtain the solution of Eq. (ITSl) in the simple 
form: 

L^Btc(£) 


da(Ea,e) = 


/c=l,2 


£ — Ea -f El E At 


( 21 ) 


The solution of Eq. m is derived from that written 
in Eq. dSU by the formal substitution a b. 

Substitution of the solutions da{Ea,e) and 
dh{Eh,e) into Eq. dH]) and subsequent integration over 
energies Ea and Eh provides the equation 

[£l - A] c(£) = ic(0). (22) 

In Eq. d211)i matrix A is defined as 
A = A-tTT ^ BjB]. (23) 

j=a,b 

The solution of Eq. (12^ can be conveniently 
written in terms of eigenvalues A^ and eigenvectors 
and of the matrix A decomposed as A = PAP^^, 


c(£) =i'^Pk 




[Pfc^c(O)] ^ 

e - Afc 


(24) 


Whereas eigenvectors pfc form columns of the matrix 
P, eigenvectors p^^ occur as rows of the inverse matrix 

p-v 

Substituting Eq. (IMl) into Eq. (1^ . coefficients 
da{Ea,e) are obtained in the form: 


2 4 


da{Ea,e)=^l^l^ 


mBip,] [p-ic( 


-■(25) 


. , . , £ — Ea + El E X^, £ — Afc 

j — 1 k—1 J 

It can be shown by direct calculation that the 
integration over Ea in Eq. m using the solution 
d2SD gives zero. Similarly, the solution for dh{Eh,£) 
obtained from the symmetry a b results in zero after 
the integration in Eq. (HSl). This confirms consistency 
of the approximation leading to Eqs. dUl) and CZl- 
Finally, using Eq. (ESI) the expression dm) for 
coefficients d{Ea,Eh,e) is expressed as: 


d{Ea,Eh,£) = 


£ — Ea — Eh + 2 El 


[itL^Btpfc] [p^-ic(O)] 
£ — Ea + El E £ ~ 


' 2 4 

EE 

. fe=i 


-|- {a 6} >. (26) 


Symbol {a 4-7- 6} in Eq. (l26ll replaces the term obtained 
by the exchange of the indicated indices. 

The inverse Laplace transform of the above 
formulas then provides the coefficients in the time 
domain. Equation (IMl) attains the following form: 


i{t) = ^ pfc[p^ ^c(O)] exp(-zAfct). 


(27) 




( 20 ) 
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Decomposition of rational functions into partial 
fractions followed by the Laplace transform leaves 
formula (ESI) as follows: 


2 4 


da(i?a,i)=EE 

j=l k=l 


[LjBtpfc] [Pfc^c(O)] 

Ea — Ei^ ~ ~ 


X [exp[-i(£’a - El - - exp(-iAfct)] . 


(28) 


Similarly, equation (1261) takes the following form in the 
time domain: 


2 4 

d{E,,Ei„t) = J2T.^ 

i=i fc=i 


[I^L^Btp,] [p^^c(O)] 
Ea + Ef) — 2 El — Afc 


+ 


exp(-iAfc<) - exp[-i(i?a - El - A5)t] 
Ea — El ~ ~ -Afc 

exp[-i(ii;a + Eb- 2EL)t] 

Eb - El + X] 
exp[-i(A;o - El- X’j)t] 

Eb - El + AJ 


+ {a -o- 6} >.(29) 


As eigenvalues A“ and Aj (A^) have positive 
(negative) imaginary parts, the above formulas 
considerably simplify in the long-time limit t —> c». 


c“(t) = [0,0,0,0]^ 

df(£;„t) = [0,0]^, j=a,b, 


d°°{Ea,Eb,t)=Y,J2 


[Pfc^c(O)] 


j=l k=l 


Ea + Eb — 2 El — Afc 


Eb — El + A^ Ea — El “f A“ J 

X ey.^[-i{Ea + Eb - 2EL)t]. 


(30) 


Both atoms are thus completely ionized in the long¬ 
time limit. In this limit, also the norm of the long-time 
joint photoelectron ionization spectrum \d°°[Ea,Eb)\^ 
is determined analytically, 


dEa / dEb\d^{Ea,Eb)\'^ = 4^' E E 
[ [liL5Btp,]*[ltL5,Btp,,] 


_ b-^j' 

Af - A‘, 


-|- {o 6} 


[Pfc^c(0)]*[p/c(0)] 


K - Afc' 


(31) 


As both atoms interact through the dipole-dipole 
interaction, energies Ea and Eb of the ionized electrons 
are correlated to some extent. This correlation is 
quantified by covariance C defined as: 


(EgEb) 

vwmy 


(32) 


In Eq. (15^ . the mean values {E^E\^) for k,l € N are 
defined as: 

(E^aEi) = I dEa I dEbE^,Ei\d°^{Ea,Eb)\^ (33) 

using the long-time joint photoelectron ionization 
spectrum \d°^{Ea, Eb)\'^ determined in Eq. ([30|). As 
these correlations have quantum origin, they can 
alternatively be described by quadratic negativity as 
it is done in the following section. 


3. Quadratic negativity for the composite 
system 


We first consider a simplified model analyzed in mi 
H5] . In this model, atom a has only two discrete 
states |0)a and |l)a whereas atom b contains the whole 
continuum \Eb) of states with energies Eb- A general 
state \tj}q) of this system described by coefficients 
dQ{Eb) and di{Eb) takes the following form: 


1^,) = J dEb [doiEb)\0)a + diiEb)\l)a] \Eb). (34) 


Quadratic negativity Nq quantifies the amount of non¬ 
separability of a state [50] as visible in negativity of the 
corresponding partially transposed statistical operator 
mii^. It is determined by negative eigenvalues of the 
partially transposed statistical operator. Quadratic 
negativity of the state \tpq) written in Eq. (l34)) has been 
derived in [48] in the form: 



-d*{Eb)dL{Eb)dl{E{,)d,{E{,)]. (35) 

Formula (l35l) can be recast as follows: 


Nq = 2j dEbQb{Eb) J dElgb{El,)n<^iEb,E{,), (36) 

where gb gives the density of states \Eb) in the 
continuum of atom b: 

1 

gb{Eb)=J2\dj{Eb)\^. (37) 

i=o 


Joint density of quadratic negativity nq introduced in 
Eq. (I36p is obtained in the form: 


nq{Eb,E'i,) = 1 - 


Elk=od*{Eb)dk{Eb)dl{E',)dq{E{,) 

gb{Eb)gb{E{^) 


(38) 


The joint density of quadratic negativity nq gives 
quadratic negativity of the bipartite system formed 
by the states {|0)a, |l)o} and {\Eb), |E^)}. According 
to Eq. (I55|) . quadratic negativity Nq is given as a 
weighted sum of quadratic negativities between the 
two-level atom a and all possible pairs of states inside 
the continuum of atom b. 
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This interpretation suggests a straightforward 
generalization to the case of two continua substituting 
|0)a \Ea) and |l)o ^ \E'^) and integrating over 
frequencies Ea and E'^-. 


A„ being coefficients of the decomposition. Functions 
fn and Qn introduced in Eq. (|46ll form the dual 
orthonormal bases. Quadratic negativity N can then 
be determined according to the formula: 


N = 



dEb 


dE'b 


[d*iEa,Eb)diEa,Eb)d*iE',,El)diE',,E',) 

-d*{Ea,Eb)d{E'^,Eb)d*iE'^,E',)diEa,E',)] .(39) 


Coefficients d{Ea,Eb) determine a state \ip) common 
to both continua at atoms a and b: 


\lP)= j dEaJ dEbd{Ea,Eb)\Ea)\Eb). 


(40) 


They also give densities Qa and g describing an 
ionized electron at atom a and both ionized electrons, 
respectively: 

g{Ea,Eb) = \d{Ea,Eb)\^, 

ga{Ea) = J dEb g{Ea,Eb). (41) 

Following the same argumentation as that applied 
earlier for the simplified model considered in |48) . 
the joint density of quadratic negativity n can be 
introduced: 


'n{Ea,E'^, Eb, El) — 2 [Pab(7?a, Eb)gab{El, E^) 

- d*{Ea, Eb)d{El,Eb)d*{El, El)d{Ea, E',)] 

X {ga{Ea)ga{El)[gab{Ea,Eb) + gab{El,Eb)\ 

X [gab{Ea,El) + gab{El,El)]}-\ (42) 


Quadratic negativity N can then be expressed using 
the densities ga and g introduced in Eq. dHJ: 

N = j dEagaiEa) J dE^gaiEl) 

X j dEb[g{Ea,Eb) + g{El,Eb)] 

X I dEl[g{Ea,El) + g{El,El)] 

X n{Ea, El, Eb, El). (43) 


Formula (H5)) can be recast into a simple form 
using the reduced statistical operator pb of the 
continuum of atom b: 


N = 2 


1 - J dEb J dEl\pbiEb,El)\ 

Pb{Eb,El) = J dEad{Ea,Eb)d*{Ea,El). 


(44) 

(45) 


An analogous formula based on the reduced statistical 
operator of atom a can also be derived. 

Another expression for quadratic negativity N can 
be reached using the Schmidt decomposition of state 
IV') given in Eq. (l40l) : 


d{Ea,Eb) =Y,fn{Ea)Xngn{Eby, (46) 

n 


N = 2 




(47) 


Any separable state gives A^ = 0, whereas we have 
N = 2{Nd — E)/Nd for the maximally entangled state 
in Nd X Nd dimensions. We note that the formula in 
Eq. (ITtI) can be rewritten as N = 2(1 — l/AT) where K 
denotes the Schmidt number of independent modes. 

We note that if atoms a and b contain also discrete 
levels the consideration of amplitude spectra with 5- 
functions allows for easy incorporation of such levels 
into the above developed description (for details, see 
[48]). 

Entanglement of two ionized electrons can easily 
be modified by filtering the energies of electrons. 
Quadratic negativity Na characterizing a state with 
energies Ea of electron a in interval {E^—AE, E^+AE) 
is obtained along the formula: 

rE°+AE pE°+AE 

Na{El)=M-^ dEaPaiEa) dEl 

JEl-AE Je°-AE 

6a.{E'a) f dEb[gab{Ea, Eb) + gab{El, Eb)] 


X J dEl[gab{Ea, El) + gablE'a, El)] 
xn{Ea,El,Eb,El), (48) 

Ma= dEa dEb\d{Ea,Eb)\\ 

Je°-ae j 

If also energies Eb of electron b are limited to 
interval (E^ — AE, E° + AE), quadratic negativity Nab 
of the resultant state is determined as follows: 

fE°+AE pE°+AE 


Nab{ElEl)=N-^ 

rEl+AE 


E«-AE 


dEa Pa {Ea) 


dEl 


lEO-AE 


Qa{El] 


Ie°-ae 


dEb[gab{Ea,Eb) + gab{El, Eb)] 


L 


eI+ae 


E°-AE 


dEl[gab{Ea, El) + PablE'a, El)] 


X n{Ea, El,Eb, El), 


(49) 


Af = 


rE°+AE 


fE'^+AE 


dEa 


lE°-AE 


Ie°-ae 


dEb\diEa,Eb)\y 


4. Long-time photoelectron ionization spectra 

In the following, we restrict ourselves to the most 
common case of two identical atoms. The discussion 
of photoelectron ionization spectra is divided into 
two parts according to the relative strength of direct 






and indirect ionization paths. For the discussion, 
we introduce a useful parametrization that generalizes 
the parametrization introduced by Fano [141 154] . In 
our parametrization, parameters /i^, flj, Vj, and Jj 
{j = a, b) are replaced by the following parameters [46] : 

7j=7r|V,P, 7j=7r|Jj-|2, 

_ - _ f^b - _ f^a 


Tj = Jj+Jj, 


+ IjQj 

^3 ~ -p 5 

^ 3 


D,j = + i)fljaL, j = a, b. (50) 


In Eq. (I50L 7 j ( 7 ^ ) gives damping of the continuum 
at atom j caused by the Coulomb (dipole-dipole) 
interaction and Fj is the overall damping. Parameter 
qj (qj) gives the ratio of indirect and direct ionization 
strengths related to the Coulomb (dipole-dipole) 
interaction. Parameter flj is linearly proportional to 
the pumping strength of ionization at atom j. As both 
atoms interact with the same laser field of amplitude 
a^, the introduction of two additional parameters is 
convenient: 


m = 


Ma’ 


n = 


2 


(51) 


Whereas parameter m gives the ratio of dipole 
moments for the auto-ionizing states at both atoms, 
parameter fl describes average pumping strength. 
A suitable parametrization is based upon common 
parameters m and 12 and parameters 7 ^, 7 ^ and qj 
of individual atoms defined in Eq. 0 . We note that 
m = 1 for two identical atoms and parameters jj, qj 
and ilj = 12 , j = a, 6 , coincide with the usual ones 
defined by Eano [Tl] . 


4-1- Spectra for comparable direct and indirect 
ionization paths 

If the values of the Fano parameters qa and qs are 
close to the unity they characterize the case in which 
both ionization paths compete each other. This 
results in the creation of the Fano zeroes at energies 
— jaqa and — ■j^qi, of electrons at atoms a 
and b, respectively. Individual spectra are formed 
by peaks that move down to lower energies with the 
increasing values of pumping parameter 12. Whereas 
the peaks are above the Fano energy for lower values 
of pumping parameter 12, they occur below the Fano 
energy for greater values of parameter 12 [14]. As 
shown in Fig. [2](a) for ionization spectrum E of 
atom a, the dipole-dipole interaction with continuum 
(described by parameters Ja and Jb) moves the peaks 
towards energies E^ and E^ of the auto-ionizing levels. 
This behavior reflects the fact that the dipole-dipole 
interaction with continuum itself ionizes atoms a and 
b and thus diminishes the contributions of direct 



(a) (Ee-E,“)/re 



(b) (Ea-E,°)/r, 

Figure 2. Long-time intensity photoelectron ionization 
spectrum la of atom a [la{Ea) = J dEb\d°°(Ea, Eb)\'^] as a 
function of normalized energy (Ea — E^)/Ta for (a) 7 a = 7 b = 0 
(independent atoms, curve without symbols), 0.1 (curve with 
*), 1 (curve with o), Jah = 0 and (b) Jab = 0 (independent 
atoms, curve without symbols), 0.56 (curve with *), 1 (curve 
with o), 1.68 (curve with A), 7a = 76 = 0; Ea = Eb = Ej^ = 1, 
IJa = Ijb = 1, 7a = 76 = 1, = 1, "I = 1- 


and indirect ionization paths whose interference is 
responsible for energy shifts of the peaks. On the 
other hand, the dipole-dipole interaction between the 
discrete auto-ionizing levels (described by parameter 
Jab) tends to form a peak close to energy (or 
E^b) in individual ionization spectra, as documented 
in Fig.[2Db). Two peaks coexist together in the spectra 
for intermediate values of parameter Jab- If the dipole- 
dipole interaction between discrete auto-ionizing levels 
is sufficiently strong, there occurs one large peak in 
the photoelectron ionization spectrum of atom a 
(or atom b). It is separated by two minima from the 
tails built by the direct and indirect ionization paths. 
The presence of two minima in the spectrum reflects 
nontrivial mutual influence between the dipole-dipole 
interaction on one side and the direct and indirect 
ionization paths on the other side. 

The dipole-dipole interaction mutually correlates 
energies of the ionized electrons at atoms a and b (see 
Fig. [3|). This results in pure states of two ionized 
electrons entangled in energies Ea and Eb- Correlations 
increase with the increasing strength of both dipole- 
dipole interactions with the continua and between the 
discrete auto-ionizing levels. In fact, there occurs 
’anti-correlation’ between energies Ea and Eb of the 
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(c) 


(a) 



(b) 


(Ea-E,°)/r, 


Figure 3. Topo graphs of long-time joint photoelectron 
ionization spectra I plotted as a function of normalized energies 
{Ea—E^)/Ta and {E^—E^)/ri) of atoms a and b, respectively, for 
(a) 7 a = 76 = 0 ? Jab = 0 (independent atoms), (b) % = 7 ^ = 1 , 
= 0 and (c) % = 7 ^ = 0, Jab = 1.68. Values of the other 
parameters are written in the caption to Fig. log stands for 
decimal logarithm. 


electrons that expresses the conservation law of energy 
in the stationary system {Ea+Ef, Ri 2El). Comparison 
of the joint photoelectron ionization spectra / plotted 
in Figs. Eljb) and (c) reveals that the stronger 
the dipole-dipole interaction the tighter the energy 
correlations. 

If the dipole-dipole interactions with the contin¬ 
uum and between the discrete auto-ionizing levels have 
comparable strengths, a well-formed minimum occurs 
in the photoelectron ionization spectrum la {h) of 
atom a (6) (see Fig. 2]). If the original spectra with¬ 
out the dipole-dipole interactions lie above the Fano 
zeros, intensity of this minimum is close to zero [see 
Fig. SKa)]. If the original spectra are below the Fano 
zeros, the minimum intensity does not reach zero but 


Figure 4. Long-time photoelectron ionization spectrum la of 
atom a for (a) Jab = 0-66 (curve without symbols), 2 (curve 
with *), 4 (curve with o), ^a = 7b = 1? ^ = 1 (b) 

7 ^ = 7 ^ = 0, Jab = 0 (independent atoms, curve without 
symbols), % = % = 1, J^t, = 0 (curve with *), % = % = 1, 
Jab = 2 (curve with o), Q = 5. Values of the other parameters 
are written in the caption to Fig. Curves with % = 7 ;, = 1 
and Jab = 2 fulfil the ’balance condition’ given in Eq. IIB.lOt . 


the intensity profile exhibits a sharp minimum accom¬ 
panied by a neighbor sharp peak [see Fig. SKb)]. Both 
cases give a clear evidence about the strong mutual in¬ 
fluence of ionization channels based on the continuum 
and the discrete levels. Conditions for the balance of 
two ionization channels based on the dipole-dipole in¬ 
teraction have been found in Appendix B using the 
Fano diagonalization approach. 

4-2. Spectra for the dominating indirect ionization 
path 

High values of the Fano parameters Pa and qt occur 
in this region, in which the ionized states are reached 
nearly exclusively from discrete auto-ionizing levels 
with energies and E^. The individual photoelectron 
ionization spectra of both atoms thus consist of 
peaks centered around energy {E^) for atom a 
(b). There occurs the Autler-Townes splitting of 
these peaks [58] for intense pumping. It occurs 
whenever the speed of populating an auto-ionizing 
level is faster than its depletion due to the Coulomb 
configuration interaction {\pa<y-L\ > |14|)- As shown 
in Fig. |S] comparing the cases of independent and 
interacting atoms, the dipole-dipole interaction with 
the continuum weakens the Autler-Townes splitting. 
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Figure 5. Long-time photoelectron ionization spectrum la of 
atom a for 70 , = 7 b = 0, Jab = 0 (independent atoms, curve 
without symbols), 7^1 = 75 = 1, = 0 (curve with *) and 

% = 7b = 0, Jab = 0.56 (curve with o); Ea = Eb = El = 1, 
qa = qb = 100, 7 a = 76 = 1, = 3, m = 1. 

This is a consequence of the fact that the dipole- 
dipole interaction makes the transfer of electrons 
from discrete excited auto-ionizing levels into their 
continua faster and so partly suppresses the effect 
of ’population reversion’ at the discrete auto-ionizing 
levels responsible for the splitting. On the other hand, 
the dipole-dipole interaction between the discrete auto- 
ionizing levels only slightly modifies the ionization 
spectra, as documented in Fig. [S] This follows from 
the considered symmetric configuration of atoms a and 
b that minimizes the influence of parameter Jab to the 
populations of discrete levels and E^. 

The dipole-dipole interaction with the continuum 
strongly modifies the joint photoelectron ionization 
spectra / that belong to the states entangled in 
energies. Dramatic change of the joint photoelectron 
ionization spectrum I caused by this interaction 
is illustrated in Fig. [6] in which the dipole-dipole 
interaction nearly completely suppresses the Autler- 
Townes splitting. 


10 log(1+l) 

I.On 


0 



(b) 

Figure 6. Topo graphs of long-time joint photoelectron 
ionization spectra 1 for (a) 7 ^ = 7 {, = 0 (independent atoms) 
and (b) 7 ^ = 7 b = 1 ; Jab = 0 and values of the other parameters 
are written in the caption to Fig. [5] 



(a) (b) 


5. Entanglement in long-time photoelectron 
ionization spectra 

The joint photoelectron ionization spectra I plotted in 
Figs. [3ljb) and (c) andlUJb) reflect strong correlations 
in energies Ea and Eb of the ionized electrons created 
due to the dipole-dipole interaction. These correlations 
emerging during the quantum evolution are non- 
classical. They reflect the bipartite entanglement that 
can be quantified by quadratic negativity N introduced 
in Sec. 3. Alternatively, they can be described by 
the von Neumann entropy S of the reduced statistical 
operators of individual electrons a and b. We note 
that entropy S' is a monotonous function of negativity 
N. Contrary to entropy, quadratic negativity N can be 
expressed via its joint density n introduced in Eq. dnj. 
This represents an important advantage as it allows to 
connect entanglement with spectral properties of the 


Figure 7. Quadratic negativity N (curve without symbols) and 
covariance C (curve with *) as they depend on (a) parameter 
7 = 7 a = 76 and (b) parameter Jab- Values of the other 
parameters are written in the caption to Fig. [2] 

ionized electrons. 

Assuming other parameters fixed, the increasing 
values of parameters Ja (7) ^^nd Jb (jb) describing the 
dipole-dipole interactions with the continua result in 
greater values of negativity N [see Fig.jTJa)]. Similarly, 
the greater the values of parameter Jab characterizing 
the dipole-dipole interaction between the discrete auto- 
ionizing levels, the greater the values of negativity 
N [see Fig. m.h)]- In Fig. m the comparison of 
curves giving negativity N and covariance C of energies 
reveals that even the classical covariance C can be used 
as a good indicator of mutual coupling of two ionized 
electrons. 
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Figure 8 . Topo graph of quadratic negativity TV as it depends 
on dipole-dipole interaction parameters 7 = 7 a =75 and J^h' 
Values of the parameters are written in the caption to Fig. [2] 


Figure 10. Quadratic negativity Na as a function of normalized 
energy {Ea — E^)/Ta of atom a; /S.E = 0.01, Jah = 1-68 and 
values of the other parameters are written in the caption to 
Fig. [21 N = 1.79. 


10 ^ 




(a) 


filtered, highly entangled states are obtained in the 
central part of the spectrum (see Fig. [10] for negativity 
No). We note that the values of negativity Nab as 
well as negativity Na depend on the length NE of the 
interval of measured energies [see Eqs. (i48l) and ((49|)]. 
The wider the interval NE the greater the values of 
negativities. However, this dependence is weak. 

6. Ionization in molecular condensates 


10’ N 

ab 



(b) 

Figure 9. Topo graphs of quadratic negativity N^b depending 
on normalized energies (Ea — ® 2 )/ra and (Eb — E^)/Tb for 
(a) 7 a =76 = 1, Jab = 0> = 0.005 [N = 0.98] and (b) 

la = ^b = 0, Jab = 1.68, AB = 0.01 [N = 1.79]. Values of the 
other parameters are written in the caption to Fig. [2] 

Whereas the dipole-dipole interactions with the 
continua and between the discrete auto-ionizing levels 
influence the ionization spectra in general in a complex 
way, they support each other in creating entanglement, 
as documented in Fig. [8] 

Spectral distribution of entanglement can be 
conveniently visualized using negativity Nab defined 
in Eq. ((49|). It characterizes a common state of both 
ionized electrons obtained by filtering the energies 
of electrons. The comparison of graphs shown in 
Figs. [9Ka) and [IKb) and also in Figs. [Hb) and [Sjc) 
reveals that the negativity N is concentrated in the 
areas with higher intensities. 

Also when the energy of only one electron is 


The dipole-dipole interaction is usually much weaker 
than the dipole interaction of atoms or molecules with 
external coherent fields. It is also much weaker than 
the Coulomb configuration interaction. This means 
that only weak modihcations of the ionization spectra 
discussed in Sec. 4 are expected in real systems. For 
example, energy shifts of the ionization peaks in the 
spectra plotted in Fig. [2] are comparable to the value 
of the dipole-dipole interaction energy. On the other 
hand, strongly correlated pairs of ionized electrons can 
be obtained in real systems provided that the process of 
ionization is sufficiently slow compared to the timescale 
characterizing the dipole-dipole interaction. 

Molecular condensates [SS] represent a typical 
example. In molecular crystals, energies of discrete 
excited states are in 1 eV, Coulomb configuration 
interaction energies in 10 meV and dipole-dipole 
interaction energies in 0.1 1 meV. Typical 

values of dipole moments are expressed in 1 D. 
Electric-field amplitudes in 1 x 10® V/m are thus 
necessary to arrive at comparable dipole and Coulomb 
configuration interaction energies needed for the 
observation of the Autler-Townes splitting. If the 
Coulomb configuration interaction energy (Va) and the 
direct dipole interaction energy (p-aCtL) have equal 
strengths, also the ionization process considerably 
slows down. Electrons then tend to occupy their 
ground states and discrete auto-ionizing states, i.e. 
states mutually interacting through the dipole-dipole 
interaction. Alternatively, the electrons can undergo 
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(a) (b) 

Figure 11. Quadratic negativity N (curve without symbols) 
and covariance C (curve with *) as they depend on dipole 
moments ji = fia = fib foi* (a) Ja = Jb = 0.001, Jab = 0 and (b) 
Jab = 0.001, Ja = Jb = 0; Ea = Eb = El = 1, fla = fj^b = 0.001, 
Va=Vb = 0.01, CL = 1. 



Figure 12. Quadratic negativity N (curve without symbols) 
and covariance C (curve with *) depending on pump-field 
amplitude cl] Ea = Eb = El = 1, /ia = = 0.001, 

Va=Vi = 0.01, ila= fLb = 0-02, Jat = 2 X 10-"^, Ja = Jb= 0. 

ionization caused by the dipole-dipole interaction 
with the continua. In both cases, the influence 
of dipole-dipole interaction to the ionization process 
dramatically increases. This results in greater values 
of negativities N [see Fig. lTlT al for the interaction with 
the continua and Fig. ITTT bl for the interaction between 
the auto-ionizing levels]. 

Considering fixed values of parameters V and fl, 
the conditions maximizing negativity N and discussed 
above can be reached choosing suitably the pump-field 
amplitude a^. This is possible because negativity N 
as a function of pump-held amplitude ol exhibits a 
well-formed maximum (see Fig. 

7. Conclusions 

Analytical solution of the model of two auto¬ 
ionization systems interacting through the dipole- 
dipole interaction has been found using the Laplace 
transform of the dynamical equations. Quadratic 
negativity together with its spectral density have 
been dehned to quantify entanglement between two 
ionized electrons developed due to the dipole-dipole 
interaction. The dipole-dipole interaction with the 


continua weakens the interference of direct and indirect 
ionization paths. This results in spectral shifts of the 
peaks in the long-time photoelectron ionization spectra 
towards the energy of the auto-ionizing state. On 
the other hand, the dipole-dipole interaction between 
the auto-ionizing states leads to the occurrence of 
an additional peak in the long-time photoelectron 
ionization spectra centered at the energy of the auto- 
ionizing state. When the auto-ionizing states play a 
dominant role in ionization and a strong pump field 
induces the Autler-Townes splitting of the ionization 
peak, the dipole-dipole interaction with the continua 
weakens this splitting. Balance of the dipole-dipole 
interactions with the continua and between the auto- 
ionizing states leads to the occurrence of a minimum in 
the long-time photoelectron ionization spectrum of an 
individual atom that resembles the Fano zero. The 
dipole-dipole interactions are responsible for strong 
anti-correlations in energies of the ionized electrons. 
They originate in the energy conservation. Both 
dipole-dipole interactions participate together in the 
creation of entanglement. The stronger the dipole- 
dipole interactions are the more entangled the state 
is. Optimal conditions for entanglement creation occur 
when the direct and indirect ionization paths are 
balanced and slow ionization is observed. ’Distribution 
of entanglement’ in the state of two ionized electrons as 
described by the density of quadratic negativity reflects 
the shape of the corresponding joint photoelectron 
ionization spectrum. 

Acknowledgments 

J.P. and A.L. gratefully acknowledge the support by 
the project LO1305 of the Ministry of Education, 
Youth and Sports of the Czech Republic. J.P. 
and W.L. also thank Operational Program Education 
for Competitiveness - European Social Fund project 
CZ.1.07/2.3.00/20.0058 of the Ministry of Education, 
Youth and Sports of the Czech Republic. A.L. 
thanks the project IGA No. PrF/2014/005. W.L. 
acknowledges Vietnam Ministry of Education Grant 
No. B2014-42-29 for support. 

Appendix A. Solution for two independent 
auto-ionization systems 

If atoms a and b are independent the quantum state 
|^/)(t) in Eq. ([5|) can be written as 

l^r"W = l^a)(t)|V’.)(t). (A.l) 

In Eq. (lA.lIl states j = a,b, describe electrons 

at individual atoms: 

= c^W|0)i+ci(t)|l)j + f dEjdj{Ej,t)\Ej). 

(A.2) 










13 


The corresponding differential equations for the 
coefficients in Eq. (|A.2I) are derived as follows: 


. d 

i — 

dt 




dj {Ej , t) 



K, 


I 

Ej - El 


Cj{t) 

10 

(A.3) 


using vectors Cj, cj = {c’q,c[) for j = a, b, and matrices 
Kj and Ij defined in Eq. ©• 

The Laplace transform leaves equations (IA.3I) in 
the form: 


[el - Kj] Cj(e) - J dEj ljdj{Ej,e) = fcj(O), 

(e - Ej + El) dj{Ej,e) - ljcj{e) =0, j = a,b. 

(A.4) 

Applying the approach similar to that used in Sec. 2 
the solution of Eqs. (|A.4I) is revealed in the form: 


cj{e) = t 

/c=l,2 


Lic,(0) 

^ + ’ 


i^j J 


i 

6 — Ej + E]^ 


E 


l]Lic,(0) 

^ + K 


j = a, b. 
(A.5) 


Matrices and eigenvalues used in Eq. (jA.5|) are 
given in Eqs. m and (EOl), respectively. The inverse 
Laplace transform provides this solution in the time 
domain: 


E I^fcCj(0)exp(fA^T), 
k^l,2 

_ llUcAO) r 

dj{Ej,t)= - ^Lxp[-i{Ea- EL)t] 

k=i.2 -El + Xi >- 

— exp(fA'it)|, j = a^b. (A-6) 


Considering the common state of both 

atoms defined in Eq. (Al), the Laplace transform of 
coefficients c{t) defined by direct product Ca{t) ^Cb{t) 
is expressed as: 


c(e)=z ^ 


[L°cj0)] [L^c,(0)] 

^ + K + K' 


(A.7) 


Similarly, the Laplace transform of coefficients 
da{Ea,t) determined by direct product Ch{t)^da{Ea,t) 
takes the form: 


da{Ea,e) = i Y 

fc,fc' = l,2 


[ltL°c,(0)] [L^Cfc(O)] 
Ea — El + A? 


£ — Ea + El + A^, £ + A“ + A^., 


(A.8) 


'k ^ ''k' J 

Using solutions (IA.7I1 and (IA.8I1 . the first term in 
Eq. (ini) can be treated as follows: 


[(£ ~ Ea + El)X + Lh] da(£'(j, £) — 


fe=l,2 


{^ — Ea + El + A^)L^da(Ea, £) — 


E 


[ltL“c„(0)][L^c,(0)] 


^ + K + >i' 


(A.9) 


k,k' = l,2 

The following relations have been used in Eq. (IA.9|) 

I„ 0 




fefe'-L'fej 


B„ = 


0 


(A.IO) 


The equality expressed in Eq. (IA.9I) implies that the 
integral term in Eq. (USD is zero when independent 
atoms a and b are considered. The same holds also for 
the integral term in Eq. o owing to the symmetry 
a O 6. We note that this conclusion can be achieved 
also by direct integration of the solution written in 
Eq. (ESI). 


Appendix B. Competition of the discrete and 
continnum dipole-dipole interactions 


Spectra in auto-ionization systems are formed by 
mutual interference of two types of ionization paths. 
One of them is based upon direct ionization originating 
from the ground state. The other uses an excited auto- 
ionizing discrete level that mediates auto-ionization of 
the system. These two ionization paths may interfere 
destructively. The presence of the Fano zero in the 
ionization spectra |14) represents a clear evidence of 
such quantum interference. Similarly, the dipole- 
dipole interaction of an auto-ionization system with 
its neighbor can be divided into two parts [IS]- One 
part influences the discrete auto-ionizing level, states 
of the continuum are modified by the other. Also 
these two parts of the interaction compete. However, 
this competition occurs at the level of population 
of the involved quantum states as the dipole-dipole 
interaction means energy transfer. This competition 
naturally modifies photoelectron spectra. The greatest 
changes in photoelectron spectra are observed when 
both parts of the dipole-dipole interaction have 
comparable strengths. Suitable conditions for this 
case can be revealed as follows. First, we apply the 
usual Fano unitary transformations of excited/ionized 
states at atoms a and b. The appropriate values of 
parameters are then found from the requirement that 
the interaction Hamiltonian iLtrans given in Eq. o 
is zero. In this case, two parts of the dipole-dipole 
interaction related to the discrete and continuum states 
of both atoms are equally strong but with the opposed 
signs. 

Invoking the Fano unitary transformations |141I54) 
on both atoms, Hamiltonians Hj written in Eq. © are 
transformed into their diagonal forms, 



(B.l) 


dEjEj\Ej){Ej\, j = a,b. 
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In Eq. (EH), states \Ej) arise from the Fano 
diagonalization, 


\E,) = f,{E,)\l), + / dE'g,iE„E')\E'), j = a,b 


(B.2) 


and 

fAE,) 


v;{EA 


9j (Ej , Ej ) 


vAEAMEj) 

Ej — Ej + ie 


+ ^{Ej 


E'). 


(B.3) 


Damping constants 7 j, 7 j = 7r|V, p, and shifted 
energies E^ E° = E°+ V J dE\Vj{E)\^/{E° -E), have 
been introduced in Eq. (IB.31) . Symbol e denotes a small 
positive number and limit e —?► 0 is assumed at the end 
of calculations. 

After the transformations, Hamiltonians Hj' 
defined in Eq. m attain the form: 




dEjg,j{Ej)aL exp(—j(0|+H.c.(B.4) 


Dipole moments p,j{Ej) describing excitation/ionization 
of states \Ej) inside the structured continua are derived 
in the form 

- / ^ — r. ^j(Ej) + Qj ^ 

d‘j{Ejj — jdj (rp \ _ ■ I (^-5) 

^3\Ej) * 

ej{E) = [E — Ej)/"fj, Qj = fXj/{irfljV*), and i stands 
for the imaginary unit. If the dipole-dipole interaction 
between two atoms is neglected, we have p,a{E^) = 
db{Eb) = 0 for E^ - E° = -"faQa and E[ - E^ = 
—"fbdb- Thus, there occurs one Fano zero with energy 
E^ in the photoelectron ionization spectrum of atom 
a and one Fano zero with energy E^ in the spectrum 
of atom b. 

These Fano zeros are concealed by the dipole- 
dipole interaction. The greatest competition of two 
parts of the dipole-dipole interaction is observed pro¬ 
vided that the matrix element of transfer Hamilto¬ 
nian iitrans between states \E^) and \E^) equals zero. 
Hamiltonian iJtrans given in Eq. o is expressed in 
these bases as 

Htrans = J dEa J dEbJ{Ea, Eb)\Ei,){Ea\ + H-C., (B.6) 

where 


J{Ea,Eb) = Jabfa{Ea)f;{Eb) 


+ J dE',j:g,{Ea,E',)fUEb) 

+ J dE',Jbfa{Ea)gUEb,E',). (B.7) 

The use of expressions written in Eqs. (IB.31) transforms 
equation (IB.71) into the form: 


J{Ea, Eb) 


JgVb 

lb[f-b{Eb) - i] 


JbVg* 

JglegiEa) + i] 


_l_ JgbVgVb - iJlVblg + iJbVglb 

lalb[^g{Ea) + i][eb{Eb) - l] 


The requirement J {E ^, E^ 


= 0 results in the 
following condition for the values of dipole-dipole 
interaction constants Jab, Jg and Jb'- 
J*gd*g , Jb^ib 


— Jab- 


(B.9) 


Ma db 

Assuming identical atoms a and b and real interaction 
constants, the condition in Eq. (IB.91) simplifies 

= -■ (B.IO) 

^ab/^ 
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